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Abstract 

We study bti production via subprocesses initiated by 6-quarks at the Large 
Hadron Collider. Both QCD and electroweak interactions are included in the 
elementary scattering amplitudes for bg — > bti. Since the additional jet in the 
final state (arising from the bottom quark) is in most cases at very low trans- 
verse momentum and very high pseudorapidity, it tends to escape detection. 
Therefore, such a process can act as a background to double-top as well as 
to single-top channels exploited in top quark phenomenology at present and 
future hadron colliders. Furthermore, if the additional 6-jet can be tagged then 
bti samples can be exploited in constraining possible effects of New Physics. 
The relevance of this reaction in such contexts is discussed and various total 
and differential rates of phenomenological interest are given. 
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1. Introduction 



The main production channel of top quarks at hadron colliders is that proceeding through 
qq and gg partonic scatterings that is, via 

qq^tt, (1) 

and 

99 -> it, (2) 

mediated by QCD, being the electroweak (EW) contributions qq — > j*,Z — > tt much 
smalleif]. 

This mode has been exploited in the recent discovery of the top quark at the Tevatron 
||. The latest value of the top mass as measured from the the CDF and DO experiments 
is m t = 175.6 ± 5.5 GeV 0, which is in very good agreement with that derived within 
the Standard Model (SM) from EW measurements at LEP, SLC, SPS, Tevatron and in 
neutrino scattering experiments [§J: m t = 177 ± 7l 19 GeV. This clearly represents a great 
success of the theory. Furthermore, it is worth recalling that when the mass of the 
boson will be determined to a precision of 50 MeV or less (at LEP 2), the EW theory will 
also significantly reduce the allowed range for the mass of the Higgs boson. 

It is clear that the discovery of the top quark and the measurement of its mass are 
only the first steps. In the near future, various details of the production process as well 
as of the decay channels will have to be thoroughly studied and other fundamental top 
parameters (the width, the couplings, etc.) will have to be determined. In this respect, a 
second mode that will be studied at the Tevatron is the production channel involving one 
top quark only ||. Although the first mechanism (i.e., double-top production) is dominant 
both at the Tevatron and the Large Hadron Collider (LHC) 0, the second (i.e., single-top 
production) is very well suited for studying the coupling between the W ± boson and the 
top and bottom quarks [[!]]. In fact, the measurement of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element Vtb cannot be easily performed via it-production and decay. On 
the one hand, top-antitop pairs are produced through a 'flavour blind' gti vertex. On the 
other hand, the dependence on Vtb entering in the width of the top is 'canceled' by the fact 
that branching fraction for the decay t — > bW^ is close to one in the SM, which predicts 
0.9989 < \V t b\ < 0.9993 at the 90% confidence level [§]. In the end, in tt samples, one 
can establish a lower limit on \V t b\ by studying the polar angle of the lepton produced 
in the top decay channel t — > bW^ — > biv^. In contrast, the single-top production rates 
are directly proportional to | V^| 2 , so that a simple measurement of total cross section can 
lead to a high precision determination of the CKM matrix element describing the mixing 
between top and down-type quarks. 

The variety of mechanisms which produce only one top in the final state is very rich. 
They all involve EW interactions. They can be conveniently grouped according to the 
final state produced and to the virtual particle content (see Ref. 0), as 

tb — production pp,pp^tbX (via s-channel W /± *-bosons), (3) 

tq — production PPiPP - * tqX (via t, w-channel H /± *-bosons), (4) 

2 They amount to « 8% of the total cross section both at the Tevatron and the Large Hadron Collider 
(see later on). 



tW ± — production pp,pp — > tW ± X (via s[t, u] -channel b[t}* quarks), 

(5) 

where q represents a light (anti) quark, b and t a bottom and a top (anti) quark, respectively. 

Although the event rate of single-top production at the Tevatron is at present rather 
smallQ, this mode will soon be investigated intensively. In fact, by 1999 the new Main 
Injector at the Tevatron should be operating, boosting the centre-of-mass (CM) energy 
\fs p p of the colliding proton- antiproton beams from 1.8 TeV of the so-called 'Run 1' (during 
the past years 1992-1996) to 2 TeV for the forthcoming 'Run 2' (1999-2001). An increase 
of the total peak luminosity is also foreseen, by one order of magnitude, up to L = 
2 x 10 32 cm -2 s _1 . A final stage of the Tevatron has also been proposed: the 'Run 3', 
always at \fs pp - = 2 TeV but with the increased luminosity L = 10 33 cm~ 2 s _1 (thus 
denominated 'TeV33', planned for the years 2003-2006). Eventually, the pp accelerator 
LHC at CERN will take over, starting operations in 2005, with a CM energy of 14 TeV 
and a peak luminosity in the range 10 32 — 10 33 cm~ 2 s _1 (Run 1). 

The processes ([[])-(§) represent the most frequent production modes of top quarks at 
the hadronic colliders of the present and future generation and they have been extensively 
studied in the literature. Many of their higher order corrections have also been calculated 



(mainly from QCD) and incorporated in the experimental Monte Carlo programs |l(J. It 
should however be noticed that the accuracy of the experimental measurements is con- 
tinuously increasing and the huge amount of data which will be collected over the next 
decade will seriously challenge the current reliability of the perturbative predictions, so 
that further theoretical efforts will soon be needed in order to match the advances on the 
experimental side. In this respect, it is clear that the increase in luminosity and in CM 
energy available will for example mean that rarer production channels of top quarks could 
yield detectable rates at the new generation of hadronic machines. 

It is the purpose of this letter to study the production of top-antitop pairs via 6-quark 
initiated processes, through the two-to-three body reaction 

bg -> btt, (6) 

at both the Tevatron and the LHC. The Feynman diagrams describing this process are 
given in Fig. la-b. In our opinion, there are at least three good reasons to expect that 
these events can eventually become interesting in top physics. 

1. Although the final state of reaction (f|) always involves an additional jet with respect 
to that produced via processes (|l|)-(||), this tends to be at very low transverse mo- 
mentum and very large pseudorapidity. Therefore, in most cases, it escapes detection 
and the final signature is precisely that of top-antitop production and decay (JTj) — (0) - 

2. As the typical (partonic) energy available at supercolliders increases, the content of 
6-quarks inside the colliding protons is very much enhanced |TT[, this yielding a much 



higher probability for 6-quark initiated scatterings to take place. 

Contrary to the case of the tt channel, for which EW contributions are due to only 
two additional diagrams, for reaction (|6|) the number of EW diagrams is much larger 
than that of the QCD ones (compare Fig. la against Fig. lb), this rendering the EW 
contributions important in the total cross section. 



3 For a review of typical single-top cross sections, see, e.g., Refs. ||, ||. 



The impact of even small corrections to ti production and decay via qq- (dominant 
at the Tevatron) and gg- (dominant at the LHC) scatterings could be crucial both in 
studying SM effects and from the point of view of the searches for New Physics. In the 
first case, bti events are interesting as they can represent a background to radiative top 
decays JT^], such as ti(g) — > bbW + W~g (which constitute a testing ground of QCD), as 



well as to three-body top decays, such as t — > bW ± Z and especially t — > bW^H |T3| 



m 



which Z, H — > bb). In the second case, one should recall that, being the mass of the top of 
the order of the EW symmetry breaking scale (the vacuum expectation value is in fact 246 
GeV) , top phenomenology represents one of the most promising places where to search for 
phenomena beyond the realm of the SM ||14]| . Therefore, to assess exactly the amount of 



all sizable corrections within the ordinary dynamics is a benchmark task to be achieved in 
order to establish possible deviations from the standard theory. 

Among the various possible effects of New Physics which could be investigated in (both 
single- and double-) top studies at hadron colliders [I5|-[]2"T|] , we discuss here only two, 



which are relevant to the case of bti phenomenology. First, various models for dynamical 
E W symmetry breaking |T5[ predict the existence of heavy colour- singlet and colour-octet 
vector states strongly coupled to the top quark, the so-called 'colorons' [|16j. These could 
produce a resonance decaying into tt-pairs. Therefore, the effects of the new particles 
should be visible in the spectrum of the ti invariant mass. Second, if there are more 
than six quarks, the CKM matrix element \Vtb\ could be anywhere between (almost) zero 
and unity, depending on the amount of mixing between the third and fourth generation, 
whereas non-universal top couplings differing from the SM form V — A should appear as 
a deviation from the branching ratio of top quarks into 6iy ± -pairs with the gauge bosons 
polarised longitudinally (which happens around 69% of the times in the SM). By studying 
the angular distributions of the charged lepton in tt-decays [BIJ , one can extract limits on 
the form of the W^tb coupling. 

A very last motivation that we put forward to justify our analysis is based on the 
observation that the three-vertex process (|]) possibly involves large (compared to the 
QCD ones, see point 3. above) contributions with two EW couplings. In some instances 
these are both W ± tb vertices (graphs 1, 5, 9 and 13 in Fig. lb), so that \V t b\ A terms enter 
in the total cross section of reaction (|6]). Indeed, we will show that diagram 5 in Fig. lb 
is the dominant EW contribution, as all the other diagrams in Fig. lb amount to only 
~ 2% of the EW part0. Therefore, provided the EW contribution is large enough in the 
total cross section (and the QCD one is subtracted), rates from process (f|) will show a 
high sensitivity to variations of \V t b\- More in general, all the couplings of the EW model 
involving a top-quark appear in the diagrams of Fig. lb. Therefore, a large variety of 
possible deviations from the SM dynamics could be tested in bti events: in particular, the 
existence of a fourth generation of quarks («4, d^) with a down-like quark with mass 
around 175 GeV. In fact, the new particle would reveal itself 'directly' in the 76^4 and 
Zd^d^ vertices of the theory. More important, when a d^ flavour is produced instead of 
a top via process (^), diagrams 1, 5, 9 and 13 of Fig. lb do not appear any longer, so 
that the total cross section should suffer from a large depletion. In contrast, in single-top 
studies via measurements of the \Vtb\ matrix element the effects of new flavours could pass 

4 Note that, for reason of gauge invariance, also diagram 1 must be included together with diagram 5 
in the computation, although this is very much suppressed because of the g — > ti splitting. In addition, 
graphs 9 and 13 are negligible, because of the splitting into ife-pairs of a real W . 



unobserved, in case of an extremely small mixing of the third and fourth generation. It 
will be one of our main concerns to quantify the relevance of the diagrams in Fig. lb 
with respect to those in Fig. la, which only involve gluon-fermion-antifermion vertices, 
so that their contribution is insensitive to the flavour. This can also be affirmed for pair 
production of ^-quarks from qq and gg-fusion, the counterpart of processes 

In summary, the final aim of our study is to assess the relevance of bti events as 
background in both double- and single-top phenomenology as well as their own importance 
in 'confirming the identity' of the particle recently discovered at Fermilab. The material 
we will present has been organised as follows. In the next Section we give some details of 
the calculation and list the values adopted for the various parameters. Section 3 is devoted 
to a discussion of the results. The conclusions are in Section 4. 



2. Calculation 

The tree-level Feynman diagrams that one needs for computing process @ are given in 
Fig. la-b. The pure QCD graphs are displayed in Fig. la whereas those involving also EW 
vertices are given in Fig. lb. To calculate the corresponding amplitudes squared we have 



used MadGraph |23[ and HELAS p4f . The codes produced have been carefully checked 
for gauge and BRS [^] invariance. The integrations over the appropriate phase spaces 
have been performed using VEGAS fl26|p|. 

The 6-quark in the initial state of reaction @ has been treated as a constituent of 
the proton with the appropriate momentum fraction distribution fb/ P (x,Q 2 ), as given by 
the partonic structure functions. So has been done for the gluon. However, as the parton 
distribution functions (PDFs) of 6-quarks inside the proton suffer from potentially large 
(theoretical) uncertainties^ (see, e.g., Ref. [|11||) and those of the gluon are not so well 
known at small x, we have produced our results in the case of several different recent next- 
to-leading structure functions, such as the sets MRRS(1,2,3) |31| and CTEQ4(HQ) [Q, 
which give excellent fits to a wide range of deep inelastic scattering data (including the 
measurements from the HERA ep collider) and to data on other hard scattering processes. 

In each case the appropriate value of A^^ (in the modified Minimal-Subtraction scheme) 
has been used. The QCD strong coupling a s entering explicitly in the production cross 
sections and implicitly in the parton distributions has been evaluated at two-loop order, 



with A-^p^ 4 calculated according to the prescriptions in Ref. and (in general) at the 



scale /i = (i.e., the CM energy at partonic level). However, since the choice of the scale 
for the structure functions and for ot s represents a source of uncertainty, we have studied 
the yields of process @ also in case of other representative values of \i. 

In the numerical calculations presented in the next Section we have adopted the 
following values for the electromagnetic coupling constant and the weak mixing angle: 
a em = 1/128 and sin 2 6y/ = 0.2320. For the gauge boson masses and widths we have taken 

5 Note that in computing the rates for process (||) we also have included the contribution due to the 
charged conjugated diagrams of Fig. la-b (i.e., those proceeding via bg- fusion). 

6 In fact, the &-sea distributions are not measured by experiment, rather these are obtained from the 
gluon distributions splitting into bb pairs by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution 
equations | j27| and, in general, the implementation of such dynamics is very different from set to set in 
those currently available on the market p8|-|34[|. 



M z = 91.19 GeV, T z = 2.5 GeV, M w ± = 80.23 GeV and T w ± = 2.08 GeV, while for the 
fermion masses we have used, in general, m& = 4.3 GeV (to match the value used in the 
MRRS(1,2,3) PDFs) and m t = 175 GeV 0. We have changed the value of m b into 5 GeV 
when using the CTEQ4(HQ) PDFs and adopted the additional values 165, 170, 180 and 
185 GeV for mt when studying the effects of a possible fourth generation of quarks. For 
simplicity, we have set the CKM matrix element of the top-bottom coupling equal to one. 

The Higgs boson of the SM enters directly in the diagrams of Fig. lb. As default value 
for its mass we have used Mh = 150 GeV, according to the best \ 2 fit as obtained from 
the analysis of the LEP and SLC high precision EW data: i.e., M H = U9±lf GeV f§. 
However, since the constraints on the Higgs mass are rather weak (a lower bound of 66 
GeV from direct searches |36j and a 95% confidence level upper limit of 550 GeV from the 
mentioned data exist [[37]) we have studied the M H dependence of the EW contributions 
of process @. 

Finally, as total CM energies of the colliding beams at the Tevatron and the LHC we 
have adopted the values y/s pp = 2 and y/s pp = 14 TeV, respectively. 

3. Results 

As a first result we quote the cross section for events of the type (H) at \fs p p = 2 TeV. 
This is very very small, around 2-3 fb. Thus, it is negligible both in double-top and 
single-top phenomenology. In fact, at lowest order, the total cross section of the former is 
around 8 pb |J whereas that of the latter is approximately 3 pb [0, |[|. That is, process (|5|) 
represents a correction of the order of one part in ten thousands. In fact, rates become even 
smaller when acceptance cuts (in transverse momentum, pseudorapidity and separation 
of the detectable particles) are implemented, in all possible top-antitop decay channels: 
hadronic, semi-hadronic(leptonic) and purely leptonic as well. The number of produced 
events is in itself tiny too, a few tens at the most at only during Run 3. Clearly, this is of 
no experimental relevance. In contrast, at the LHC, the total cross section is more than 
three orders of magnitude larger. Therefore, in the following we will focus our attention 
to the case of the CERN hadron collider only. 

Before presenting our results for the LHC a few words are needed concerning the 
dynamics of process (|6]). In particular, we would like to point out that its cross section is 
finite over all the available phase space (when all masses are retained in the calculation). 
This allow us to safely calculate the rates of process @ in both the following cases: (i) 
when the additional 6-quark in the final state is produced in the detector acceptance region 
and (ii) when it escapes detection either because at low transverse momentum or because 
at large pseudorapidity. Indeed, from the point of view of top-antitop phenomenology the 
latter case is more interesting. In fact, as mentioned in the Introduction, a reason to study 
reaction ([]) is that it represents an irreducible contribution to top-antitop production and 
decay when the additional jet goes undetected. The former case is instead interesting as 
it can represent a background to radiative top decays and to three-body top decays. 

Furthermore, we should also mention that in calculating the differential distributions 
involving the top decay products we have in general resorted to a Narrow Width Approx- 



imation (NWA), as we have written the heavy quark propagator as 



- m t / p \ a 
p 2 — in 2 + im t r \T t ) ' ^ 

where T t is the tree-level top width and where the numerical value adopted for T has been 
10 -6 . This way, one is able to correctly reproduce the rates for ttX production times the 
(squared) branching ratio [BR(t — ► 6T4^ ± )] 2 , as it is done in many of the experimental 
simulations, which indeed do not include finite width effects [j3gfl . Incidentally, this sim- 
plifies considerably the numerical calculations as one can integrated out the Breit-Wigner 
dependence of the two top quarks, thus reducing by two the number of variables of the 
multidimensional integrations. (We have verified that the spectra of quantities involving 
the 6's, the leptons and the light jet produced in the top decays suffer little from NWA 
effects.) The only exception has been made when plotting the distribution in the partonic 
CM energy, for which the exact results have been obtained (i.e., T = T t ). 

Finally, we identify the jets in the final state with the partons from which they originate, 
we introduce no jet energy smearing and apply all cuts at parton level. 

3.1 Theoretical error 

The rates of process (|6]) depend strongly on the 6-quark parton distribution function inside 
the proton. On the one hand, such density is very poorly constrained by the data, as even 
the extreme case in which this is neglected the derived parton densities (such as the GRV 



PDFs |39], f5J) can adequately fit the available data. On the other hand, there has been 
a lot of theoretical activity in the recent years in the field of heavy quark distributions 
(although mainly focused to the case of charm quarks) |3T|, [34], f|l], In the latter 

context, it should be mentioned that two approaches had dominated the literature in 
the past years: the so-called TFN-scheme, where TFN stands for Three-Flavour-Number 
(such as in the GRV sets), and the FFN-scheme, where FFN stands for Four(and Five)- 
Flavour-Number (such as in the MRS and CTEQ parton densities). In the first case, 
only the 'light' flavours {g,u,d,s} participate in the parton dynamics inside the proton, 
and heavy flavours can only be created in scattering processes (i.e., 'flavour creation'). In 
the second case, also the 'heavy' quarks {c, b, (t)} are numbered among the initial state 
partons, provided the energy scale of the interaction is large enough (though in most cases 
the 'massive' parton are indeed treated as 'massless', once the density is turned on at the 
appropriate threshold^) , so that heavy quarks in the final state can also be obtained by 
exciting the corresponding flavour inside the proton (i.e., 'flavour excitation'). It turns 
out that (see, e.g., Refs. |j42| ) the TFN-scheme is the most suitable for the heavy quark 
component of the PDFs near threshold whereas well above this regime is the FFN-scheme 
that should be used. 

More recently, a consistent formulation of heavy flavour dynamics within the 
perturbative-QCD (pQCD) framework has been given, in Refs. [42| . The new 

treatment encompasses both the 'flavour creation' and 'flavour excitation' mechanisms 
and is valid from the heavy quark threshold up to the high energy regime. It reduces 
to the two above approaches in the appropriate limits: i.e., to the TFN scenario when 



7 For a comparison between mass-dependent and mass-independent evolution of PDFs, see Ref. |0 



/i ~ m q and to the FFN one if \x 3> m q (where q = c, Although the three formulations 
in Refs. PH , |4l| , 42| | are slightly different, their approach is basically the same. Further- 



more, based on the improved theory for heavy quark dynamics, global analyses and new 
PDF packages have been made available, such as the MRRS(1,2,3) |3l| and CTEQ4(HQ) 
|J§ sets. 

Preliminary comparisons between the various 'heavy quark sets' have been performed 
in Refs. [[y], fH, |43|, though results are not conclusive yet, since similarly fitted parton 
distributions in the various schemes are not available in the literature at the moment, to 
allow for a consistent study |34| . Being such a comparison beyond the scope of this paper, 



we confine ourselves to the computation of some relevant rates of process (|6|) with the 
four new sets of parton distributions. The total cross section for bg — * bti, before any 
acceptance cut, is given in Tab. I for MRRS(1,2,3) and CTEQ4(HQ), at the LHC. Indeed, 
Tab. I shows that the differences among the four sets of PDFs are reasonably contained. 
In fact, if one assumes as default the MRRS(3) value, the largest difference occurs with 
respect to CTEQ4(HQ) one, about +10%. Rather than an absolute value, such a difference 
should be considered as a lower limit (for the present time) on the error affecting the rates 
of process (|6|) due to the PDFs. 

Before proceeding further, it should be noticed that, together with process ©, one 
should also consider another reaction, namely 

gg — > bbti (via g — > bb splitting), (8) 

which gives the same final state as the process bg — > bti if one of the 6-jets is close enough 
to the beam pipe. As a matter of fact, the two scattering processes @ and (||) (in which 
one of the two incoming gluons splits into bb pairs) should be considered simultaneously 
and their cross sections summed up with a subtraction of a common part in order to avoid 



double counting (see, e.g., Ref. |34| , [411). That is, from the sum of the 'flavour excitation' 
scattering (||) and the gluon-fusion mechanism (|8|) , one has to subtract the piece due to the 
configurations in which the quark in the internal 6-line (from the g — *> bb splitting) in (§) is 
on-shell and collinear to the incoming gluon. One would expect that the total results for 
the combination of bg — > bti and gg — > bbti events should be less //-scale dependent then 
the rates for process (||) alone JOJ. As a criteria to decide whether the additional 6-jet 
falls inside the detectable region (thus mimicking a two-to-three body hard scattering) or 
not, we use the cut prip) < 20 GeV on either of the bottom quarks (see eq. flT0D below). 

The total rates due to processes (||) and (|J) with the mentioned subtraction as a 
function of the //-scale are reported in the left-hand column of Table. II. From there, one 
can argue that the error associated with the scale dependence of reaction ([|) should be 
below 30% or so. In fact, this is the difference between the value of the cross section, say, 
at \i = 2m t and at 1 TeV. Such numbers correspond to the case of the MRRS(3) PDFs 
(our default in the following), though we have verified that similar rates also occur for the 
other three sets adopted here. 

Finally, for completeness, we also have calculated the contribution from events of the 
type 

gg — > bbtt (via g — > tt and g — > gg splittings), (9) 



8 As this new approach effectively interpolates between the preceeding two, it is often referred to as the 
Variable-Flavour-Numbcr (VFN) scheme. 



in which one of the gluons splits into either ti or gg pairs, with one of the final state 6's 
missed along the beam line. As a matter of fact, also such events contribute to the total 
rate for bti production, on the same footing as the @ ones do. The corresponding cross 
sections (with one of the fe's having pr{b) < 20 GeV) as a function of the /i-scale can be 
found in the right-hand column of Tab. II (i.e., in brackets). They amount to roughly 18% 
of the total rates in the nearby column. 

For sake of simplicity in the numerical calculations, in the remainder of the paper we 
will only consider the rates of the process bg — > bti, thus ignoring for the time being its 
interplay with those of reactions @ and (H). 

3.2 bg — > bti phenomenology 

Contrary to the case of the Tevatron, at the LHC effects due to process (H) can be per- 
ceptible in various instances. For starting, the cross section is rather large per se, as it 
amounts to 8 pb. Therefore, with 10 fb _1 of Run 1 at the LHC, it yields some 80,000 bti 
events per year. In one selects semi-(hadronic)leptonic SM decays, then the event rate is 
around 24,000. 

It is interesting to separate the total cross section given previously into its pure QCD 
component (i.e., diagrams in Fig. la) and that involving EW interactions as well (i.e., 
diagrams in Fig. lb). This is done in Fig. 2, for the LHC. In addition, to establish whether 
the presence of the Higgs graphs 4, 8, 12 and 16 (in Fig. lb) has any influence on the total 
cross section we vary Mh in the range between 60 and 500 GeV. As one might expect such 
contributions to be particularly relevant near the H — > ti threshold (see diagrams 4 and 
16 in Fig. lb), we have enlarged the interval 340 GeV ~ Mh ~ 360 GeV in the insert in the 
middle of the plot. As a matter of fact, the effect is visible in the EW contribution but 
imperceptible in the total cross section. Thus, Higgs diagrams are of no special concern, 
whichever the actual value of Mh is. This implies that no further theoretical uncertainty 
related to the unknown parameter Mh enters in our results. 

The important feature in Fig. 2 is that the EW contribution is rather large in the total 
results. At the LHC and for Mh = 150 GeV this is around 3 pb, yielding some 30,000 
events during Run 1 (which become 9,000 after implementing the semi-hadronic(leptonic) 
SM branching ratios of the two top quarks). This is around 40% of the total rate of 
process (|J). If a fourth generation down- type flavour gLl of mass m^ 4 ps 175 GeV is instead 
produced, these rates change dramatically. As previously mentioned, this is mainly due to 
the absence of the dominant EW diagrams 1 and 5 in Fig. lb. However, also the different 
EW couplings between neutral gauge bosons (i.e., 7 and Z) are responsible for differences. 
Tab. Ill shows the total cross section (pure QCD graphs plus the EW ones) for events of 
the type bg — ► bd^d^ against that of bti, for five values of (conservatively compatible 
with the latest measurements of m t ). One can see that the rates are very different. Indeed, 
those for bd±d± production practically coincide with the QCD contribution of process @ 
only (compare the rates in Tab. Ill with those in Fig. 2, when m t = 175 GeV). This 
pattern remains unchanged for all masses considered. 

Clearly, if one can isolate top-antitop events compatible with an additional 6-quark in 
the final state, then it could be possible to search for deviations such as those revealed 
in Tab. III. This could be achieved in the semi-hadronic(leptonic) top decay channel by 
requiring two 6-tags and one high p? lepton and by reconstructing the resonant fermion 



mass via a three-jet system. The requirement of two vertex tags imposes an additional 
reduction of the detected rates. By assuming an efficiency of 60% per fiducial 6-jet |44|, [15 



the overall one to tag two of these would then be 36%. Of the original bti events, some 
8,000+3,000 would survive during Run 1, further reduced by a factor of 9 after the following 
detector acceptances (on two 6-jets)0: 

|t?(&)|<2, p T (6)>20GeV, 

\r}(£)\ < 2.5, pr(£) > 10(20) GeV, 

\ V (j)\ < 4, p T (j)>20GeV, 

> 20 GeV. (10) 



T 

The final numbers are approximately 880 and 330 for the final rates from the diagrams in 
Fig. la-b, respectively. This means that roughly 1200 events are expect within the SM 
from bti production in the detection channel bbijjj, where £ — e, \i and j is a jet. Because 
of the dominance of diagrams involving two W^tb couplings in the E W component of bti 
events, several hundred of bit events will show a \Vtb\ A dependence, which could possibly 
be exploited in experimental analyses to furnish an independent new measurement of the 
CKM matrix element. A fourth generation d-type quark with mass similar to that of 
the top would only produce the equivalent of the QCD component of bti, that is, 38% less 
events. Finally, we have verified that no intrinsic difference exists in the kinematics of the 
QCD and EW components in the spectra of the lepton and jets which could be profitably 
exploited to separate these two parts of the cross section. 



3.3 Double- and single-top phenomenology 

The numbers in Tab. I and Fig. 2 certainly compare rather poorly to the total rates of 
both double- and single-top events. In fact, the cross section for process (||) amounts to 
1% of the former and 2% of the latter. 

In the first case, this clearly implies that effects due to events (^) can hardly be dis- 
entangled in ti phenomenology. Indeed, one should recall that the total cross section for 
pair production of top quarks is made up by the two components (0)-(|) only, each of 
these yielding the same signature and each of these with rates much larger than those of 
©. Of the 800 pb at the LHC, 730 pb come from gg- and 70 from gg-fusion ||. How- 
ever, if particular kinematic configurations are selected to reduce the leading order ti rates 
(e.g., to disentangle some New Physics or when studying the mentioned radiative and/or 
three-body top decays), then bti events could become important. 

Fig. 3 (left-hand side) plots the invariant mass distributions of the system tiX for 
events of the type (|I])-(0) and (^|). From there it is clear that if one looks for inclusive 
signals (i.e., X represents any additional particle in the final state) and studies the total 
invariant mass of the final state, then a 'resonance' (dashed curve) should appear some 
180 GeV above the ti threshold at 2m t . Furthermore, note that also the peak in the 
M t i spectrum produced by events of the type @ (dotted line) is shifted by several GeV 
(upwards) with respect to that due to qq- and gg-fusion (solid line). 

9 In parentheses is the pt threshold for charged leptons used as trigger. 



The cosine of the polar angle of the (anti)lepton from the two top decays (signal of a 
possible contamination of V + A coupling in the W ± tb vertex) in btt samples is significantly 
different from that originated by processes (0)-(@), as can be appreciated in right-hand 
side of Fig. 3. As the coefficient of a possible V + A chirality violating current could well 
be at the level of few percent compared to that of the V — A term, it is important to 
recognise and subtract the contribution of btt events from the experimental sample. 

If one compares the rates due to process © with the yield of each of the subprocesses 
contributing to single-top production, one discovers that the 400 pb cross section (including 
all charged conjugated processes) for pp — > tX (X does not include here a second t-quark) 
is built up in the following way ||: pp — > tqX yields 253 pb, pp — > tW ± X produces 128 
pb whereas pp — > tbX gives 19 pb ||. Therefore, the rates of process (|6]) are similar to 
those for single-top production in association with a bottom quark. 

The latter process has a particular relevance in single-top phenomenology. In fact, 
there are three parton subprocesses yielding final states of the type, say, tbX (i.e., for 
top-antibottom pairs). Namely, (i) q'q — > tb; (ii) q'g — > tbq; (iii) q'q — > tbg (see Ref. 
Among these, the 'ideal' one to be exploited in phenomenological studies would be q'q — > tb 



46] . This is because the corresponding cross section can be reliably calculated, for several 
reasons. First, the light (anti) quark distribution functions inside the proton are evaluated 
at moderate values of x where they are well known. Second, the QCD corrections to this 
process can be calculated easily (they plug into a pure EW process) up to the order 0(a 2 a ). 
Those from the initial state are taken into account by simply constraining the quark- 
antiquark flux via a measurement of the rates for q'q — > Iv^. Those from the final state 
can be incorporated without ambiguities as there are no collinear and soft singularities. As 
a matter of fact, the two radiations do not interfere until the 0(a 2 s ) order. Furthermore, 
the kinematics is that of a two-to-two body process, whereas cases (ii) and (iii) involve an 
additional light jet in the final state, thus case (i) is easier to reconstruct experimentally. 
More in general, the process q'q — > tb can boast the same advantages with respect to any 
of the other single-top production channels. For one or more of the following reasons. 
Firstly, the latter always involve either a gluon, a b or both in the initial statefj^. Secondly, 
they proceed via QCD interactions. Thirdly, singularities can occur at higher orders in a s . 
Therefore, in the first case, they suffer from larger uncertainties due the rather unknown 
corresponding PDFs and, in the second and third cases, the inclusion of higher order QCD 
effects is in general much less trivial^. 

The cross section for q'q — > tb at the LHC, before any acceptance cuts, is around 9 pb 
[[46] (including C.C. processes), almost the same as that of process (Q), thus the latter must 
be considered as possible dangerous source of background events. To establish its relevance 
in this context, we proceed as follows. First, we calculate the cross section of events of the 
type (§) in which the additional 6-quark in the final state escapes detection. In this case, 
the final state recorded in the detectors is ti. Then, we select semi-leptonic(hadronic) 

10 Note that in case (ii) only graphs with gluons coupled to light quarks are considered, as those involving 
g — > bb and g — > ti splitting form a separate gauge invariant set and need to be considered apart, because 
of the definition of the 6-sea quarks in the partonic distributions (see discussion in Ref. ||). 

n The only exception is qq — > tW~b, which is however two orders of magnitude smaller H and with a 
more complicated final state. 

12 This is also true for the simple two-to-two EW process q'b — > qt, where collinear and infrared diver- 
gences along the massless q'q line spoil the advantages from the factorisation of the strong corrections at 
the order 0(a s ). 



decays tt — > bbW + W~ — > bbiuejj, where £ = e,// and j represents a light quark jet, 
originated by the fragmentation of a u-, d-, s- or c-quark. To appreciate the kinematics 
of the additional final state 6-quark in reaction (|J) one can refer to Fig. 5, where its 
distributions in transverse momentum pt and pseudorapidity 77 are plotted, along with 
those of all the other final state particles, as they would appear at the LHCQ. We introduce 
here a distinction between the three 6's in the final state of (|6|), that we will maintain in 
the following. This because they have different kinematic behaviours. We call 'prompt V 
that produced in association with the tt pair, 'direct V that coming from the decay of the 
top in Fig. la-b and 'indirect V the remaining one from the antitop^. Such a distinction 
has phenomenological relevance for the 'prompt V for the key role that it plays in the 
determination of the observed rates of process @, whereas in the other two cases it has 
been made for sake of illustration only. In the first case, one can appreciate the striking 
behaviour of the parton in Fig. 5, as this tends to be emitted almost exclusively along the 
beam pipe, with px < 20 GeV and \rj\ > 4. For the other two b's, one should recognise 
that, although the production dynamics of the 'direct V and the 'indirect V is different 
(at least in the diagrams of Fig. lb, not in those of Fig. la) and this can be spotted in the 
left hand side plots of Fig. 5 (and 4 too), the differences cannot possibly be disentangled 
experimentally. We shall now proceed by performing the same selection procedure outlined 
in Ref. j|6| and will compare the yield of process (^) with those presented there. 



The acceptance cuts implemented to simulate the detectors are those already mentioned 
in eq. (]T0|), supplemented by the further requirements of separation (in general, on the 
'direct' and 'indirect' b's): 

\AR bb \ > 0.7, \AR M \ > 0.7. (11) 

Note that in applying cuts (|10|) (|TT|) to process (^j in the context of single-top phenomenol- 
ogy one only has four jets in the selected final state (two fe-jets and two light flavour ones), 
as one of the three produced 6-quarks (in general, the 'prompt' b) is assumed to be lost 
along the beam pipe. The cross section of bg — > bti events after the above acceptance 
cuts is approximately 1600 fb. After implementing the branching ratio 2/9 for electronic 
and/or muonic decays of the tagged W ± one gets o « 360 fb, thus 3,600 events during 
Run 1 at the LHC (with ten inverse femtobarns). The rate obtained for process q'q — > tb 
in Ref. |46j after the same requirements and in the same detection channel is (reading 



from Tab. 1 in that paper) 580 fb (including charge conjugation). The signal is thus a 
factor of approximately 1.6 above the background from process ([]). The application of 
the cut M b i > 110 GeV (see again Ref. fl46"f ) imposes a factor of two of reduction on bti 
events. Finally, we obtain an additional drastic rejection against the latter, because of the 
additional W ± in the final state. In fact, as can be appreciated in Fig. 5, any additional 

13 Note that we also have plotted (in Fig. 4) the same rates for the Tevatron, although we have already 
assessed that they are of no relevance at such a collider. We have done so simply to show quantitatively 
that the percentage of events (^ with the first 6-quark produced escaping detection because of the finite 
coverage in pr and 77 is much higher at larger energies: compare the solid curves in the two plots on the 
left hand side of Figs. 4 and 5. In contrast, the detector requirements on the other particles in the final 
state are equally effective at both Tevatron and LHC energies (dashed and dotted lines in the left windows 
and all curves in the right ones). 

14 The labelling works the other way around for the last two b's if bg induced diagrams are considered. 
Note that in the caption of Figs. 4-5 we have applied the distinction between 'direct' and 'indirect' also 
top the (anti)top decay products. 



light flavour jet j in reaction (EH falls inside the acceptance region defined by pr(j) > 20 



GeV and \f]{j)\ < 4 |46 |. We estimate that the final cross section for bti events, after light 
flavour jet rejection and after including vertex tagging efficiency would be of a few femto- 
barns, that is one order of magnitude smaller than the signal and two orders smaller than 



the background rates obtained from W ± bb, W ± jj, W^Z, tbj and tt background events 



46f . Therefore, also in single-top SM phenomenology events of the type (g) are well under 
control, though they represent a sizable correction (of several percent) to the q'q — > tb 
signal, which should be considered when searching for New Physics. 



4. Summary and conclusions 

In this paper we have studied top-antitop production in association with an additional b- 
quark at the Tevatron and the LHC The interest in this process for top studies at hadronic 
colliders comes from the fact that the 6-quark is in most cases at low transverse momentum 
and large rapidity, so it tends to escape the acceptance region of the detectors. Therefore, 
such events enter naturally in the candidate ti sample and should then be considered if 
one wants to carry out detailed studies of top quark properties. 

We have shown that at the Tevatron the total cross section for such events is negligibly 
small, around 2-3 fb for the upgraded Fermilab pp collider with \/s^ = 2 TeV. In contrast, 
at the CERN accelerator the corresponding number is 8 pb. This should yield about 80,000 
events after Run 1 for an accumulated luminosity of 10 inverse femtobarns. Though such 
a number is certainly affected by a large indetermination due to the fe-structure functions, 
we believed it reasonably large anyway so to consider in detail the relevance of btt events 
at the LHC, both as a background to double- and single-top events and as a signal on 
its own. As for the error associated to the parton density of 6-quarks, we have given an 
estimate of its lower limit, of about 10%, by comparing the total rates as obtained from 
four very recent sets of parton distributions, all implementing the heavy quark dynamics 
in the context of the newly developed Variable-Flavour-Number (factorisation) scheme. 
Furthermore, we have established a 30% uncertainty of the total rates, depending on the 
choice of the factorisation scale /i, which was varied between twice the top mass and the 
TeV scale. (Indeed, to minimise the impact of the two mentioned errors, we have neglected 
the absolute normalisation of the relevant differential rates.) 

In the case of double-top physics, it has been shown that the total rate of btt events 
is small, as it represents a correction of 1% only to the QCD signal via qq, gg — > ti. 
However, this effect could be above the experimental accuracy and further enhanced by 
dedicated selection procedures (especially in order to disentangle possible effects due to 
New Physics), so that it should be considered when proceeding to MC simulations. In 
particular, the kinematics of the top quarks in btt events is very different from that of the 
leading ti production: for example, in the Mux invariant mass, where X represents any 
additional particle in the final state. In fact, the three-jet mass distribution is much 
broader then the M# one and the peak of the resonance is shifted by +180 GeV whereas 
the two-jet one M u - is indeed similar to that of ti events but the shape of the threshold 
resonance is distorted. Finally, the spectra in the polar angle of the decay leptons differ 
significantly from those produced by ordinary ti events. 

In the case of single-top physics, it has been demonstrated that the bti production 



rates are similar to the yield of the process q'q — > tb, which has been advocated as the best 
channel to probe the W ± tb vertex of the underlying theory. Once an appropriate selection 
strategy of tb events is adopted, the bti ones are however greatly reduced. Nonetheless, the 
latter should anyway be included in the experimental simulations aiming to measure the 
Cabibbo-Kobayashi-Maskawa matrix element \V t b\ and to test the vector/axial structure 
of the top-bottom-PF 1 * 1 coupling, as violations of the SM dynamics could well be at the 
same level and show the same kinematic features as those due to the bti corrections. 

As for signal on its own, events of the type bg — > bti have been proved to be extremely 
sensitive to the possible presence of a fourth generation of quarks, involving a down-type 
fermion d<± with mass similar to that of the top, which could then mimic the latter. This 
is due to the fact that the EW contribution (dependent on the flavour of the produced 
fermion) is almost 40% of the total bti rates and almost coincides with diagrams involving 
off-shell W^- currents, which are naturally absent in case of bd±d& production. Thus, if the 
measured cross section suffered from a large depletion, this would indicate that the heavy 
particle detected recently at the Tevatron is indeed a new flavour of the theory. Conversely, 
about 40% of the total bti Standard Model cross section shows a \V tb \ A dependence, so that 
top-antitop events accompanied by an additional 6-quark could be exploited to obtain a 
new independent measurement of this crucial quantity. 
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Table Captions 

[I] Total cross sections (QCD and EW summed) for process @ at the LHC for four 
different sets of structure functions. Errors are as given by VEGAS. 

[II] Total cross sections (QCD and EW summed) for process (||) at the LHC for four 
different choices of the scale parameter. In parentheses are the corresponding rates 
for process (|9]). The MRRS(3) structure functions have been used. Errors are as 
given by VEGAS. 

[III] Total cross sections (QCD and EW summed) for processes bg — > bd\d± (see the 
text) and (|6]) at the LHC for five different masses m t /md 4 . The MRRS(3) structure 
functions have been used. Errors are as given by VEGAS. 



ure Captions 



Lowest order Feynman diagrams describing process (§): (a) the 0(a^) contribution; 
(b) the 0(a s al m ) contribution. The package MadGraph [^] has been used to pro- 
duce the PostScript file. In (b) 'A' represents a photon and the dashed line identifies 
the SM Higgs boson. 

Cross sections for process (||) at the LHC as a function of the Higgs mass. The 
MRRS(3) structure functions have been used. Solid line: process (|), 0(af) + 
0{a s a 2 em ) rates. Dashed line: process (j^j, 0(a s a 2 em ) rates. 

Left window: differential distributions in invariant mass of the systems ti in events of 
the type (solid line), bti in events of the type (|6|) (dashed line), ti in events 

of the type (|6]) (dotted line). Right window: differential distributions in polar angle 
of the lepton in events of the type (|l])-(@) (solid line), antilepton in events of the 
type (^-(0) (dashed line), lepton in events of the type (H) (dotted line), antilepton 
in events of the type @ (dash-dotted line). Normalisations are to unity. 

Differential distributions in transverse momentum (upper two plots) and in pseu- 
dorapidity (lower two plots) of the particles in the final state of process (H) at the 
Tevatron. All QCD and EW contributions have been considered. On the left hand 
side, spectra of the three 6-quarks: 'prompt V (solid line); 'direct V (dashed line) 
and 'indirect V (dotted line). On the right hand side, spectra of the leptons and light 
quark jets: 'direct leptons/jets' (solid and dotted lines) and 'indirect leptons/jets' 
(dashed and dot-dashed lines). Each curve is normalised to unity. The MRRS(3) 
structure functions have been used. 



Same as Fig. 4 at the LHC. 



crtot (fb) 


PDFs 


A* = 14 TeV 


MRRS(l) 


8117 ± 25 


MRRS(2) 


8135 ±26 


MRRS(3) 


8101 ±24 


CTEQ(4HQ) 


8996 ± 23 


M H = 


150 GeV 


no acceptance cuts 



Tab. I 



crtot (fb) 


\i (GeV) 


= 14 TeV 


2m t 


10474 ±38(1856.3 ±8.6) 


400 


10030 ±36(1719.2 ±6.9) 


500 


9459 ±32(1505.8 ±6.2) 


600 


8913 ±30(1356.4 ±5.6) 


700 


8539 ±27(1242.3 ±5.1) 


800 


8258 ±24(1153.4 ±4.9) 


900 


8019 ±22(1082.5 ±4.6) 


1000 


7802 ± 19(1010.3 ±4.4) 


M H = 150 GeV 


no acceptance cuts 


MRRS(3) 



Tab. II 



(Jtot (fb) 


m d Jm t (GeV) 


bg — > bd^di 


bg — > btt 


165 




7112 ±27 


11303 ±34 


170 




6142 ± 25 


10062 ±27 


175 




5319 ±25 


8101 ±24 


180 




4626 ± 23 


8093 ± 20 


185 




4029 ± 17 


7305 ± 18 


M H = 150 GeV 


no acceptance cuts 


MRRS(3) 






w = 14 TeV 





Tab. Ill 



Diagrams by MadGraph 




2 



Fig. 1a 



Diagrams by MadGraph 




graph 13 graph 14 graph 15 



Diagrams by MadGraph 




graph 16 



Fig. 1b 
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